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Resume  of  Recent  Research 


L.  Ridgway  Scott 

oo 

L  Estimates  for  Finite  Element  Methods 

Two  projects  have  been  completed  in  the  1980-81  academic  year.  The  first, 
a  collaboration  with  C.  Goldstein  begun  several  years  ago,  concerned  estimates 
for  various  methods  for  the  Dirichlet  problem  on  a  domain  with  a  curved 
boundary.  The  paper  developed  a  central  body  of  estimates  for  the  approximate 
Green's  function  for  arbitrary  methods  satisfying  certain  assumptions.  These 
were  then  augmented  by  special  estimates  for  three  specific  methods,  yielding 
optimal-order  error  estimates  for  each  method  in  the  maximum  norm.  A 
manuscript  entitled  "Optimal  maximum  norm  error  estimates  for  some  finite 
element  methods  for  treating  the  Dirichlet  problem"  has  been  submitted  for 
publication  to  CALCOLC . 

The  second  project  was  a  collaboration  with  R.  Rannacher  begun  in 

August,  1980.  We  were  able  to  settle,  to  some  extent,  an  issue  concerning 

the  rate  of  convergence,  in  various  norms,  of  the  finite  element  method  using 
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piecewise  linear  elements.  We  were  able  to  show  that  the  error  is  0(h) 

(h  =  mesh  size)  in  any  L  ,  1  <  p  <  »  ,  and  0(h)  in  the  Sobolev  space 
W1  for  any  p  in  the  range  1  <  p  <  °°  .  The  best  results  previously 

p  - 

had  included  a  factor  of  |log  h|  to  a  power  depending  on  p  .  For  example, 

the  best  rate  previously  known  in  w^  was  0( h  | log h | )  .  Our  research  left 

00 

open  the  rate  of  convergence  in  L  ;  an  example  indicates  that  the  known 
estimate  0(h  |loghj)  may  be  best  possible.  A  manuscript  describing  these 
results,  entitled  "Some  optimal  error  estimates  for  piecewise  linear  finite 
element  approximations"  has  been  submitted  to  Mathematics  of  Computation  and 
is  available  as  MRC  Technical  Summary  Report  #2191. 

-  1  - 


I 


Numerical  Simulation  of  Nonlinear  Dispersive  Waves 

A  numerical  method  for  solving  the  equation 
(*)  +  au  +  8u  u  -  yu  -  5  u  =0 

t  X  X  XXt  XX 
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was  developed,  analyzed  rigorously  and  used  to  compare  the  predictions  of  | 

the  model  equation  (*)  with  laboratory  experiments,  jointly  with  W.G.  Rritchard 
and  J.L.  Bona.  A  paper  describing  this  work,  entitled  MAn  evaluation  of  a 
model  equation  for  water  waves"  will  have  appeared  in  PhJA.  Trans .  Roy.  Soc. 

Lond.  A  vol.  302 ,  pp.  157-510. 

The  previously  mentioned  numerical  method  was  also  used  to  compare  solutions 
of  (*)  to  those  of  the  so-called  Korteweg-de  Vries  equation.  The  numerical 
scheme  was  used  to  integrate  (*) ,  starting  with  initial  data  for  which  the 
exact  solution  to  the  Korteweg-de  Vries  equation  was  known.  This  work 
will  be  described  in  a  manuscript  currently  in  preparation,  jointly  with 
Bona  and  Pritchard,  which  studies  more  generally  the  similarities  between 
(*)  and  the  Korteweg-de  Vries  equation. 


Numerical  Models  for  Plasticity 


A  finite  element  method  using  discontinuous  piecewise  polynomials  was 
developed  and  proved  to  be  convergent,  jointly  with  C.  Johnson.  A  report 
on  this  work,  "A  finite  element  method  for  problems  in  perfect  plasticity 
using  discontinuous  trial  functions,"  has  appeared  in  Nonlinear  Finite 
Element  Analysis  in  Structural  Mechanics ,  W.  Wunderlich  et  al^  eds. , 
Springer-Verlag,  Berlin,  1981,  pp.  307-324.  Research  concerning  rates  of 
convergence  of  the  method  has  been  done  recently  by  Johnson. 
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Numerical  Methods  for  the  Transport  Equation 


A  collaboration  with  J.  PitkHranta  is  in  progress  concerning  rates  of 
convergence  of  the  so-called  discrete  ordinates  method  for  the  transport 
equation.  The  simplest  form  of  this  occurs  in  the  so-called  "slab"  geometry, 
namely , 


(J|j)  f 

P  —  (x,p)  +  <Ji(x,p)  =  a  }  i|)(x,v}dv  +  f(x) 
i|/  ( 0 ,  p )  =  0  for  p  >  0  ,  tf>(l,U)  =  0  for  p  <  0 


Here  x  e  {0,1]  and  ;j  8  [-1,1]  .  The  discrete  ordinates  method  consists 
in  discretizing  the  integral  in  (**)  via  numerical  quadrature  and  the 
remaining  spatial  differential  operator  by  a  finite  difference  or  finite 
element  method.  A  key  feature  of  the  problem  is  the  lack  of  smoothness  of 
\Ji  near  p  =  0  .  A  manuscript  is  currently  in  preparation  that  deals  with 
the  problem  (**} ,  giving  rates  of  convergence  in  terms  of  both  the  p-mesh 
and  the  x-mesh  approximation  parameters.  In  particular,  results  are  given 
for  Gaussian  quadrature,  which  is  the  most  widely  used  method  for  the  p-integral 
discretization.  These  results  are  based  on  estimates  for  Gaussian  quadrature 
recently  obtained  jointly  with  R.A.  DeVore.  The  latter  work  is  described  in 
"Error  estimates  for  Gaussian  quadrature  and  weighted-L^  polynomial  approximation," 
which  will  have  appeared  as  an  MRC  Technical  Summary  Report  and  will  be  submitted 
to  a  journal. 


Computation  of  Fluid-Flow  Problems 


During  the  *80- '81  academic  year,  a  numerical  computer  code  has  been 
under  development  at  MRC,  jointly  with  W.O.  Pritchard  and  Y.Y.  Renardi.  One 


of  the  ultimate  goals  is  a  code  for  a  problem  with  a  free  surface.  Various 
hurdles  have  been  leapt  toward  this  goal,  but  a  working  code  for  the  free 
boundary  problem  is  not  yet  available.  Subroutines  for  calculating  the  free 
surface  from  certain  quantities  and  for  computing  the  exact  solution  to  the 
Jeffery-Hamel  flow  problem  have  been  written  and  tested.  (The  Jef fery-Hamel 
flow  is  to  be  used  as  a  test  problem;  it  has  the  property  that  the  nonlinear 
term  does  not  vanish  identically.)  Also,  a  program  for  flow  in  a  pipe 
has  been  written  and  tested.  The  latter  has  been  written  in  modular  form, 
so  the  extension  to  more  complicated  flow  regions,  including  a  free  surface, 
should  be  relatively  easy. 


Resume  of  Recent  Research 


Mitchell  B.  Luskin 

Research  has  been  undertaken  in  four  areas  during  the  past  year.  A  manu¬ 
script  has  been  prepared  with  Professor  William  Martin  and  Mr.  Leonard  Lorence  of 
the  Department  of  Nuclear  Engineering  at  the  University  of  Michigan  on  a  finite 
element  discretization  of  the  neutron  transport  equation.  The  iterative  solution 
of  the  resulting  linear  system  by  a  block  Gauss-Seidel  method  has  also  been  analyzed. 
This  procedure  has  been  shown  to  require  less  storage  and  fewer  arithmetic  operations 
than  the  direct  solution  by  Gaussian  elimination. 

A  manuscript  has  been  prepared  with  Professor  Rolf  Rannacher  of  the  Univer¬ 
sity  of  Erlangen-NUrnberg  (West  Germany)  which  demonstrates  that  for  parabolic  equa- 
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tions  with  rough  initial  data  a  0(k  /t  )  convergence  rate  for  the  Crank-Nicolson 
time  discretization  scheme  can  be  obtained  unconditionally  if  an  appropriate  smooth¬ 
ing  procedure  is  performed  (such  as  adding  A  backward  Euler  steps).  Further,  a 
report  has  been  prepared  with  Professor  Ivo  Babuska  of  the  University  of  Maryland 
on  an  adaptive  time  discretization  procedure  for  parabolic  problems  which  has  been 
published  in  the  Proceedings  of  the  Fourth  1MACS  International  Symposium  on  Computer 
Methods  for  Partial  Differential  Equations.  A  more  extensive  paper  on  this  subject 
is  currently  being  prepared. 

Finally,  a  paper  has  been  written  with  Professor  Blake  Temple  of  Rockefeller 
University  which  proves  the  existence  of  global  weak  solutions  to  an  initial-boundary 
value  problem  for  a  nonlinear  hyperbolic  system  which  models  fluid  flow  in  a  pipe. 
The  effect  of  friction  has  been  modeled  by  adding  a  quadratic  zero  order  term  to 
the  system  of  conservation  laws  for  compressible,  frictionless  flow.  A  priori 
bounds  have  been  obtained  by  means  of  a  nonincreasing  functional  that  is  compatible 
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with  the  friction  effects.  The  boundary  values  for  this  problem  cannot  be  imposed 
weakly,  so  new  results  on  the  regularity  of  the  solution  at  the  boundary  have  been 
given. 


Technical  details  of  the  research  may  be  found  in  the  following  publications 
of  Mitchell  B.  Luskin: 


On  the  smoothing  property  of  the  Galerkin  method  for  parabolic  equations  (with  R. 
Rannacher) ,  SIAM  J.  Numer .  Anal.,  to  appear. 

On  the  existence  of  global  smooth  solutions  for  a  model  equation  for  fluid  flow  in 
a  pipe,  J.  Math.  Anai.  Appl.,  to  appear. 

On  a  finite  element  method  to  solve  the  criticality  eigenvalue  problem  for  the 
transport  equation  (with  J.  Descloux) ,  SIAM  J.  Numer.  Anal.,  to  appear. 

An  adaptive  time  discretization  procedure  for  parabolic  problems  (with  I.  Babuska) , 
Proceedings  of  the  Fourth  IMACS  International  Symposium  on  Computer  Methods  for 
Partial  Differential  Equations,  1981. 

The  existence  of  a  global  weak  solution  to  the  nonlinear  waterhammer  problem  (with 
B.  Temple),  Comm.  Pure  Appl.  Math.,  submitted. 

An  analysis  of  the  Crank-Nicolson  scheme  for  parabolic  problems  (with  R.  Rannacher)  , 
in  preparation. 

An  iterative  procedure  for  the  finite  element  solution  of  the  multidimensional 
transport  equation  (with  W.  Martin  and  L.  Lorence) ,  in  preparation. 
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